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Abstract 



The analysis of tagged B s — > J/ i/j(f> decays determines the CP phase <fr s in B s — B s mixing 
i-C ' with a two-fold ambiguity. The solutions differ in the sign of cos 4> s which equals the sign of the 

width difference AT S among the two B s mass eigenstates. We point out that this ambiguity can 
be removed with the help of B s — > D S K decays. We compare untagged and tagged strategies and 
find the tagged analysis more promising. The removal of the sign ambiguity in AT S can be done 
with relatively low statistics and could therefore be a target for the early stage of B s — > D S K 
studies. 

> 

ff) \ The theoretical description of B s — B s mixing involves the elements Mf 2 and T\ 2 of the mass and 

decay matrices, respectively [1]. The precise measurement of the mass difference AM, = Mh — Ml 
between the heavy and light mass eigenstates of the B s system determines |Mf 2 | [2]. Yet theoretical 
uncertainties still permit new-physics contributions to |Mf 2 | of order 30%. The new contributions 
O ■ to Mf 2 can even exceed the Standard Model value in magnitude, because Standard-Model and new- 

physics contributions generally come with an arbitrary relative complex phase. To probe new physics 
in Mf 2 further, it is therefore mandatory to study the phase of Mf 2 experimentally. The width differ- 
ence between the heavy and light mass eigenstates of the B s system is given by 



Mf 



Ar s =r L -r H = 2|r s 12 |cos(/> s with s = arg(-^J. (1) 

Since T\ 2 is unaffected by new physics, measurements of (j) s probe new physics in Mf 2 . It is useful to 
decompose (f> 3 as <j) a = 4>^ M + 4>f with the two terms denoting the Standard-Model and new-physics 
contributions to <fi s , respectively. Currently most experimental information on <f>^ stems from the 
decay B s — > J/ip<p. The CP phase in this decay is the difference <f>f — 2(3 S between arg Mf 2 and twice 
the phase of the b -> ccs decay amplitude with (3 S = arg[-V ts V* b /(V cs V* b )} = 0.020 ± 0.005 = 
1.1° ± 0.3°. While it is safe to neglect (/>f M = (4.2 ± 1.4) • 10" 3 = 0.24° ± 0.08° [3,4] and 
to identify (f> 3 with (f)f-, we keep 2/3 s non-zero in our formulae. The untagged decay ^bI — > J/ipcf) 
provides information on AT S cos((f>f — 2(3 S ) and | sin(^ — 2/3 s )| [5]. These measurements have 
been combined with experimental constraints on the semileptonic CP asymmetry a s \ [6] to determine 
the allowed ranges for cj)f [4,7]. Recently, the CDF collaboration has presented a tagged analysis of 
B s — > J/t/jcf) [8] with the two solutions 

^-2$, € [-1.36,-0.24] or <j)f - 2/3 s G [-2.90, -1.78] @68%CL. (2) 

The quoted ranges correspond to the analysis in Ref. [8] which constrains AT S in Eq. (Q} with the 
theoretical value |T 12 | = 0.048 ± 0.018ps _1 [3,4,9]. Eq. ([T]) implies signAT s = sign cos0 S) 
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so that the two solutions in Eq. © correspond to AF S > and Ar s < 0, respectively. Lifetime 
measurements in the components of the angular distributions of an untagged B s — > J/ip4> sample 
determine |Ar s |, which implies a four-fold ambiguity in cf>f [4, 10, 11]. Neglecting the small S for 
a moment, all quantities which can be extracted from B s —> J/ tfxfi and also a s \ suffer from the same 
two-fold ambiguity 4>f <-> ir — (f>f visible in Eq. ©. Thus at present we do not have any information 
on the sign of Ar s . The two solutions in Eq. (O correspond to different values of cos^i^, where 
Si and 62 are strong phases. In order to resolve the ambiguity in sign cos (j) s = sign AT S one must 
determine sign cos (5^2. If this is done with naive factorisation [12], the second solution in Eq. © 
is obtained. However, if the strong phases measured in Bd —>■ J/ipK* [13] are used, one finds the 
first solution in Eq. © (see the discussion in [8]). It should be noted that the SU(3)f symmetry links 
Bd —> J/tpK* only partially to B s — > J/ipcj). Only the component of the (j> meson with U-spin 
equal to 1 belongs to the symmetry multiplet of the K* . The decay amplitude into the equally large 
U-spin-zero component cannot be related to Bd — > J/ipK*. Since there is also no reason to trust 
naive factorisation in B s — > J j we conclude that the sign ambiguity in AT S is unresolved. 

The tagged decays ^b\ — > DJK^ were proposed to determine the angle 7 of the unitarity triangle 
(UT). They do not involve any penguin pollution and are therefore hadronically very clean [14]. How- 
ever, these decays are also sensitive to a possible new phase in B s — B s mixing and really determine 

4>^ — 2f3 s + 7 (up to a tiny correction of order 0.1°). An exhaustive study of { bI — > DJK^ including 
the effects of a non-zero (fif can be found in Ref. [15], which also focuses on the determination of 7 
assuming that 4>f has been determined unambigously with other methods. In this paper we propose 
to view B s — > DJK^ from a different angle: We exploit that 7 is well-measured at B factories and 
show that B s DJK± can be used to discriminate between the two solutions with cos 4> s > and 
cos 4> s < in Eq. ©. This information can be found with relatively low statistics, much before studies 
of this decay mode at LHCb become competitive for the determination of 7. The resolution of the sign 
ambiguity in AT S can be achieved either with a lifetime measurement in untagged B s —>■ DJK^ 
decays or by inspecting the sign of the oscillating term in a tagged B s — > DfK^ data sample. We 
will discuss both strategies below. Since the CDF experiment has already gathered more than 100 
B s — > DJK^ events [16], the determination of sign AT S maybe even within reach of the Tevatron. 

7 is well-known from the decay B^ — > p + p~ , which measures the UT angle a plus a potential phase 
of new physics in Bd—B ( i mixing: If the measured CP asymmetries in Bd — ► p + p~ and Bd — > J/ipK s 
are combined to solve for 7 = n — a — (3, any new physics in Bd — Bd mixing drops out. Exploiting 
the smallness of the penguin pollution in Bd —>■ p + p~ and using QCD factorisation [17] Ref. [18] 
finds 7 = 71° ± 5°. For simplicity we define 

7s = 7 - 2& = 69° ±5°. (3) 

Alternatively one can include other B — > pp modes and use isospin symmetry to control the penguin 
pollution [19]. In principle this analysis comes with discrete ambiguities for 7 as well. However, the 
global analysis of the UT only permits the one solution for 7 quoted above. 

The time-dependent — > DJK ± decay rates involve [1, 14, 15] 

q (D-K+\B s ) 

Ad ^ k+ P(D7K+\B_ S ) 

q {D+K-\B s ) 
DfK~ ~ p{DtK -\ Bs) 



\^D7K+\ e 



-i( ls+< t>*-6) 



1 c -i(ls+^+S)_ 
K D7K+\ 
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Here q/p encodes B s — B s mixing in the usual way and 5 is a strong phase which equals zero if the 
matrix elements are computed in the factorisation approximation [15]. The B s — > DJK^ 1 decays 
are colour-allowed tree level processes and the factorisation approximation is exact in the limit of a 
large number N c of colours. We point out that the only 1/N C corrections to the matrix elements stem 
from annihilation topologies which are empirically known to be small. The remaining corrections to 
the large-iV c limit are quadratic in 1/N C , see e.g. [20]. Of course the full-flesh tagged analysis can 
determine 6 [14, 15], but for our purposes it is sufficient to know that 5 is small. We conservatively 
assume \5\ < 0.2. 



We introduce the shorthand notation 



2|A 



D7K+\ 



i + IW+l 2 



(4) 



and note that < b < 1. For realistic values |A 



D7K+\ 



0.4 one finds b « 0.7. 



The time-dependent decay rates for the four relevant processes are [1] 

T (B s (t) Df K±) = iVe" rfl '[cosh(^^)± (l-6|A D - x+ |)cos(AM s t) 



Ar t 

- b cos(7 s + <j£ T 5) sinri(— ^) + b sin( 7s + <j)f =f 8) sin(AM s t) 
T (B s (t) -► DjK±) = Ne- rst [cosh(^^) T (l -b\\ D7K+ \) cos(AM s t) 

at t 

- b cos(7 s + <j)f^5) sinh( — ) - b sin( 7s + <pf =F S) sin(AM s t) 



(5) 



Here T s = (Tl + Tjj)/2 and N is a normalisation constant. Now the untagged decay rate for the 
D^K T reads 



decay mode ^bI 



r [DfK*^] = r (B s (t) -> dJk±) + r (B s (t) -> d?k±) 



2Ne 



-r s t 



cosh( 



b cos( 7s + 0^ =F <5) sinh( 



(6) 



which is just the familiar two-exponential formula with the time-dependent factors exp[— T^t] and 
exp[— Tjjt]. In practice one can determine two quantities from the untagged decay rate in Eq. ©, the 
branching fraction and the lifetime measured in the considered mode ^bI — ► DfK T . The normalisa- 
tion constant N can be related to the CP-averaged branching fraction [1, 11]: 



b(b ] s 



DJK±) 



B(B S - D?K±) + B(B S - DfK±) 



ivr. 



q - (AT s )2/4 



l-6cos(7 s + <^f T<5) 



at 5 

2T, 



(7) 



From Eq. (O one finds 



sin(7 s + 0^) sin 5 Ar s 



(■ft) 



DfK-) + B(B' s ^DjK 



1 — cos (7 S + 0^) cos 5 2r s 



(8) 



Once this ratio of branching fractions is measured at the level of a few percent, it will be useful to 
place tighter bounds on | sin (7,, + cj)f) sin<5| and may help the tagged analysis. 
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Figure 1: The tagged ( bI — > J/ijxp analysis gives one solution for sin(<^ — 2j3 s ) corresponding to 
the two-fold ambiguity <f>f — 2(3 S *-* tt — 4>f- + 2(3 S . The left plot shows how the measurement of 
L in Eq. (TTOb can resolve this ambiguity For sin(</>^ — 2f3 s ) < the upper branch corresponds to 
cos 4>f~ > and AT S > 0, while the lower branch corresponds to cos < and AT S < 0. For 
sin((f)f — 2j3 s ) > the situation is reversed. The right plot shows the coefficient S of the tagged 
analysis. The upper curve is for cos 4> s > meaning AT S > 0, the lower curve corresponds to 
cos (j) s < meaning AT S < 0. Both plots are for b = 0.7 and \S\ < 0.2. The blue (dark) curves 
correspond to 7 = 71°, the red (medium grey) and green (light) curves correspond to 7 = 66° and 
7 = 76° , respectively. 



For our purposes we need the lifetime information: A maximum likelihood fit of a time evolution 
given by Eq. © to a single exponential oc exp[— ^ D f K ±t] determines [11,21]: 

AT 

r D ^± = T s + b cos( 7s + <f£ T S) = T s + b cos( 7s + <f>f T cos <f£ |Tf 2 1 , (9) 

where we neglected corrections of order (IS.T s ) 2 /T 2 S . Comparing the rates T D + K _ and T D - K + gives 
the same information as Eq. ®. More important for us is the average of the two widths: Defining the 
quantity L through 

TDtK ~ \ T ° JK+ ~ r s = 6cos,5cos(7s + 0f)cos0f|r s 12 | = L \F S 12 \ , (10) 

one first realises that the dependence of L on 5 is only quadratic, so that the uncertainty from S is 
inessential in view of the error on 7 S in Eq. ©: \5\ < 0.2 implies 0.98 < cos 5 < 1. Second we 
verify from Eq. (fTOb that we can resolve the ambiguity in by comparing the lifetime measured 

in — ► D S K with 1/T S , provided that (pf differs from or tt. This feature is illustrated in the 
left plot of Fig. [T] For example, the central values in the two intervals in Eq. (O both correspond 
to s'm((frf - 2j3 s ) = -0.72. But the solution with AT S > comes with L > 0, while Ar s < 
implies L < 0. We do not recommend to fit the data to a single exponential, because Eq. Q is only 
correct, if e.g. the experimental acceptance does not vary with the decay length. Instead we propose 
to determine with the exact formula in Eq. © with AT S expressed as AT S = 2|rf 2 | cos (f)f and 
|rf 2 | = 0.048 ± O.OlSps- 1 [4]. Further T s could be fixed to the theoretical value T s ~ l/r Bd . 

Next we discuss the tagged analysis: With T(B s (t) D S K) = [T{B s (t) -» D-R + ) + T(B s (t) -» 
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DfK )]/2 we encounter the CP asymmetry 

T(B a (t) -> D S K) - T(B s (t) -> 

T(B s (t) + T(5 s (t) D S K) 

6 cos <5 sin( 7s + <ftf) sin( AM a t) 

cosh(Ar s /2) - 6cos*cos(7 a + ^)sinh(Ar s /2) 

The coefficient of the oscillating term, 

5 = b cos S sin (7^ + (f) s ) (12) 

also permits the removal of the discrete ambiguity in 0^, because the replacement of <fif — 2f3 s by 
n — 4>f + 2(3 S changes S dramatically, as shown in the right plot of Fig. [TJ S even discriminates 
between the two cases <f>f = and <fif~ = it, which are the two possible cases in the class of new 
physics models without new sources of CP violation. 

Comparing the untagged and tagged method we find that discriminating the two branches for L in the 
left plot of Fig. [T] means a lifetime measurement with an accuracy of roughly 2% requiring at least 
2500 events, because the difference of the two solutions for L|rf 2 |/r s hardly exceeds 0.04 and even 
vanishes if is close to or n. The tagged measurement looks better, even though tagging costs 
roughly a factor of 12-20 in statistics. The two solutions with sin(<^ — 2j3 s ) = —0.72 correspond 
to S ~ 0.3 and S ~ —0.6 and a fairly small data sample should permit to discriminate between the 
two solutions. Finally we remark that one can eliminate 5 altogether, if both L and S are measured 
precisely: Eqs. (fTOl and (fl2|) combine to 

S 

tan(7 s + <f)f) = — cos^>f . (13) 

Now Eq. ([T3l has four solutions for <pf and two of them can be eliminated with the information on 
the sign of S. The remaining two solutions are not related by cf)f — 2(3 S <-> it — <j>f + 2(3 s , so that 
in combination with B s —> J/^<p the discrete ambiguity is lifted with the help of Eq. (fT3l . We 
emphasize that Eq. (fT3l has been discussed before in Ref. [15], where tan(7 s + <f)g) is expressed 
in terms of the coefficient of sinh(Ar s t/2) in Eq. ©. The extraction of this coefficient has a sign 
ambiguity if the sign of AT S is unknown. Through Eqs. (flOl and (fT3l we have merely expressed 
tan(7 s + (j)^ ) in terms of \T\ 2 1 and cos to eliminate the implicit dependence on sign AT S . 

In conclusion we have discussed the removal of the two-fold ambiguity in the extraction of from 
tagged B s —> J/ip<p decays. We have shown that B s —> DfK T data can be used to resolve this 
ambiguity. This analysis can be done with relatively low statistics, well before B s — > DfK T decays 
become competitive for the determination of 7. Comparing untagged with tagged analyses we find 
that the tagged analysis is more promising, despite of the penalty from small tagging efficiencies at 
hadron colliders. 
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